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' Abstract. - We extend the idea of quantum phase transitions of light in the photonic Bose-Hubbard 

model with interactions to two atomic species by a self-consistent mean field theory. The excitation 
of two-level atoms interacting with coherent photon fields is analyzed with a finite temperature 
C ' dependence of the order parameters. Four ground states of the system are found, including an 

isolated Mott-insulator phase and three different superfluid phases. Like two weakly coupled su- 
perconductors, our proposed dual-species lattice system shows a photonic analogue of Josephson 
effect. The dynamics of the proposed two species model provides a promising quantum simulator 
for possible quantum information processes. 
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Driven by quantum fluctuations at absolute zero temperature, quantum phase transitions 
CO ' (QPTs) have been intensively studied in interacting many-body problems [T]. Typically, it 

, is difficult to control and probe such exotic quantum phenomena for strongly correlated 

electronic systems in condensed matter physics. Optical lattices - artificial crystals made by 
interfering laser beams - offer a versatile platform for studying the QPTs of trapped Bose 
gases [51 . In this situation, one can describe the many-body dynamics from a Mott-insulator 
^> ', (MI) phase to a superfluid (SF) phase in a gas of ultracold atoms with periodic potentials by 

\^ • using the Bose-Hubbard model that includes an on-site two-atom interaction and hopping 

5^ \ between adjacent sites [3]. 

Instead of weakly interacted ultracold atoms, photons are non-interacting bosons, and in 
purely photonic systems there is no possibility to have any QPTs. For a pure Bose system, 
the conducting phase at zero temperature is presumably always superfluid [J. However, en- 
gineered composites of optical cavities, few-level atoms, and laser lights can form a strongly 
interacting many-body system to study the concepts and methods in condensed matter 
physics from the viewpoint of quantum optics. In this case, a photonic condensed matter 
analogue could be realized with state-of-the-art photonic crystals embedded with high-Q 
defect cavities. Like optical lattices for matter waves, the quantum phase transitions of 
photonic insulator (excitations localization) to superfluid (excitations delocalization) are 
predicted by the Bose-Hubbard model with additional photon-atom interactions [SHZ]- Re- 
lated quantum transitions also have been predicted in a Heisenberg spin 1/2 Hamiltonian 
[8], two species Bose-Hubbard model [9], and solved exactly in the one dimensional case 
[10) . This opened the possibilities to study critical quantum phenomena in conventional 
condensed matter systems by manipulating the interaction between photons and atoms. 
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Fig. 1: Illustration schematic for the proposed system. An array of high-Q electromagnetic cavities 
is formed in the configuration of the square lattice, and each cavity contains a single two-level atom 
of the type A or B, which is spaced at intervals. In this configuration, the center cavity has four 
nearest inter-species neighbors (different atomic type) and four next-nearest intra-species neighbors 
(the same atomic type), and eg are the transition energies for the atomic species A and B, 
respectively. The incident optical field interacting with the bipartite lattice is also shown in the red 
color. 



Recently, we illustrate the generality of the method by constructing the dressed-state 
basis for an arbitrary number of two- level atoms (TLAs) [TT]. As the number of TLAs 
increases, collective effects due to the interactions of atoms among themselves give rise to 
intriguing many-body phenomena. With the Dicke-Bose-Hubbard Hamiltonian, we show 
that the Mott insulator to superfluid QPTs with photons can be realized in an extended 
Dicke model for an arbitrary number of two-level atoms. As quantum many-body effects 
associated with Bose-Einstein condensations in optical lattices have been reported for a 
long-time [l^ , it is believed that the realization of strongly correlated many-body photonic 
cavity quantum electrodynamic(QED) systems in experiments is also likely to happen soon. 

From the view point of a practical experimental setup, instead of modifying the Q-value 
of individual optical cavities, it is more easier to allocate different atomic species in different 
photonic cavity sites. More recently, based on scanning electron microscopy technologies, the 
addressability in a two-dimensional optical lattice for ultracold atoms has been demonstrated 
experimentally down to the single-atom and single-site level |13| . A natural consequence 
of the studies on the QPT of light is to consider the atom-light interacting system with 
different atomic species. In this work, we investigate theoretically QPTs of light from the 
MI to SF quantum phase transitions in high Q-value optical cavities with two species of 
atoms based on the photonic Bose-Hubbard model embedded with the Jaynes-Cummings 
interaction Hamiltonian. A rather rich low temperature phase diagram emerges even for a 
relatively simple structure of the cavity lattice. As the phase diagram of two-component 
bosons on an optical lattice [14l[T5], we show that separated phases corresponding to MI 
and composite SF exist. Moreover, these composite SF phases have an analogue to the 
well known Josephson effect in a two weakly coupled superconductor. We expect that more 
controllable light-wave technologies should lead to an enhanced understanding of the QPTs 
of light with distinctive properties, improvements in the organization of the ground-state 
wave function, and the introduction of new applications. 

For two atomic species interacting with photons, we consider an array of two-dimensional 
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square photonic bandgap micocavities with two types of two-level atoms (TLAs) labelled 
as A and B, as illustrated in FigHJ The array is made by high-Q electromagnetic cavities. 
Each cavity contains a single TLA of the type A or B, which is spaced at intervals. The 
atoms are assumed to interact strongly with photons in the photon-blockade regime so that 
a coherent composite polaritonic (photon and atom) system is formed in the presence of 
radiation fields [5l[ll]. Photons can travel from one cavity to another. We assume that 
travelling process for photons is dominated by the hopping between adjacent sites. The 
value of the hopping matrix elements is a function of the distance between the cavities, 
which is a tunable parameter in our system. 

The interaction between a two-level atom and the quantized photon mode of an optical 
cavity is described by the Jaynes-Cummings model |16j . with the two-site Hamiltonian 
for the i-th unit cell (two cavities numbered as 2i and 2i + 1 for A- and _B-type atoms, 
respectively), 

Trtwo— site t I z , / t I t — \ 

+ '^«2i+l«2i+l + eB,2i+lO'B^2i+l + .9-B,2j:+l(a2i+lO'B_2j + l + "L+i^b.2j:+i)'(1) 

where tA{B),i is the transition energy for the TLA in i-th cavity and uj is the radiation 
field frequency which we have assumed to be the same at all cavities. The atom-photon 
coupling, gA{B),i: is assumed to be real here. a| and are the raising and lowering operators 
for photons, respectively. '^'a[b) i excites the A[ov i?)-type two-level atom from the low 
energy state to high energy state and '^a{B) i describes a reverse process. By including the 
hopping process for photons, the total Hamiltonian describing our proposed system for a 
configuration of 2N cavities is 

27V N 

^ = - E + E ["^(^^2. + n2^+l) + Hl"''^] , (2) 

where k^j is the hopping matrix element and i, j = 1, . • . , 2N are the cavity indices. We 
have also introduced a chemical potential term /i for photons which control the overall 
strength of the photon field, and Ui = a\ai is the photon number operator. For simplicity, 
in this configuration, each cavity has four nearest inter-species neighbors (different atomic 
type) with the hopping coefficient ka^^b — t^B^A = , and four next-nearest intra-species 
neighbors (the same atomic type) with the hopping coefficient ka-it^a = i^b^^^b = k' . 

To solve the model Hamiltonian, we introduce a self-consistent mean field approach by 
approximating 

A similar decomposition is also applied to the terms associated with a^cr^^^)- It is convenient 
to introduce two superfluid order parameter ipAiB) for photons and two TLA order parameter 
Ja(b) for atomic species A and B, respectively, where 

^PA^{a2^)^{al), (4) 
iPb = (a2i+i) = (4i+i)> (5) 
Ja(B) = (f^!l(s)> = (^A(B))- (6) 

The order parameter introduced here is indeterminate at the point of the phase transition 
in our photonic Bose- Hubbard Hamiltonian [SHT]. The system is in the photonic superfluid 
phase for a non-zero order parameter ipA{B) 7^ 0, while in the insulator phase for a zero order 
parameter ipA — ipB ~ 0. And Ja(b) can be viewed as an order parameter for the atomic 
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coherent states. The physical meaning for Ja{b) ss a coherent state order parameter for 
the atoms Ues on the coherent superposition of the ground and exacted state of the atoms. 
When the system temperature goes down, the coherent superposition of the ground and 
exacted state of the atoms becomes larger, from Eq. pUj) . and it reaches in its largest value 
at absolute zero temperature as expected. With these order parameters, the possibility to 
have QPT at low temperature would be demonstrated in the right parameter regimes. 

In the following we assume that both ipA{B) and Ja{B) real numbers and spatially 
independent, i.e. no spontaneous symmetry-breakings in the translation and rotation for the 
infinite system considered here. The assumption that ipAiB) a-nd Ja{b) are real numbers will 
be justified later by the self-consistent mean-field solution. In general it is very difficult to 
diagonalize the model Hamiltonian in Eq. ([2]) and find out all the desired QPTs of the sys- 
tem. However, it may be instructive to diagonalize part of the total Hamiltonian first. With 
the mean- field approximation, Eq.® can be decomposed into two Hamiltonians, labelled 
as and Rp, where the atomic and photonic operators can be considered separately. For 
one atom (A or B) inside each cavity, the mean-field Hamiltonian for atomic operators is 



H 



'a(B) = <^A(B)<y\(B)+ 9A(B)'>Pa(B) 



^A(B) + ^A(B) 



(7) 



This Hamiltonian can be diagonalized easily in the form 



H 



f 



A(B) - ^A(B) 



1 



I 



(8) 



with the eigen-energy 



Ea{B) = { \9A{B)i>A{B) 



^A{B 



1/2 



(9) 



Then the TLA order parameter J a(b) can be derived as 



Ja{i 



'^a(b)3a(b) 



IE 



A{B) 



Tanh[- 



E 



A{B) 

2T 



(10) 



where we introduce the temperature for the system, T, by setting j3 = l/ksT with the 
Boltzmann constant fc^ = 1. Note that here Ja{b) is a real number, and consistent with 
our initial assumption. From Eq. (|10p . it can be seen clearly that only at zero temperature, 
the TLA order parameter Ja{b) has a maximum value. 

For the radiation fields, the related photonic mean-field Hamiltonian in Eq.© can be 
separately for the two sites 2i and 2i + 1 as 



N N 



"-2i+l 



N 



i=l 

N 



i=l 



(11) 



'^id"-l(^3 + Y[i^- M)"2j+1 + gB.2'i.+lJB{a2i+l + 4j+l) • (12) 



These Hamiltonians can be diagonalized by the Fourier transforms for the even- and odd- 
numbered sites, i.e. 



dA.k 



a-B,k 



N 



N 



^eJ"'=-^^-+iaB,2,+i. 



(13) 
(14) 
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In terms of a^(B).fc: the combined photonic mean- field Hamiltonian of two sites in a unit 
cell becomes 



HP 



- RP + RP 



(15) 



where 



Oo(fc) = -4k' Cos{k^) Cos{ky) + uj - fj,, 
ni{k) = -2k [Cos{k^) + Cos{ky)] , 



(16) 
(17) 



with the Fourier vector k = kxX + kyij. The assumption that Ja{B) is spatially independent 
implies that the atoms couple only to the k — photon mode. In other words, we only 
search for the spatially homogeneous solutions in our mean-field theory. 

The photonic mean-field Hamiltonian in Eq. (jisp can be diagonalized by introducing 



(18) 
(19) 



symmetric and antisymmetric photon field operators a^y^.k SiJid aant,k, i-G- 

0'sym,k = {aA,k +aB,*:)/V2, 
aant,k = {aA,k — aB,k)/V2. 



And in terms of asym.k ^md Uant.k field operators, the photonic mean-field Hamiltonian HP 
becomes 



fc 

+ QsymJ sym{(^\ym,k=0 ^sy7n,k=o) + gantJant{o^ant,k=Q O,ant,k=o), 



(20) 



where 



^sym (^) 



r!o(fc) + r!i(fc), 



gsymJsym = —f^{9 aJ A + QB J b) , 



3ant J a 



V2' 
1 

V2 



{qaJa - qbJb)- 



(21) 
(22) 

(23) 
(24) 



Here the stability of photon fields requires a non-zero energy (frequency) , i.e. ^sym{ant) (k) > 
0. Equivalently we have the condition {uj — fi) — 4k' > 4k required for photon fields. 

From the coupling of photon field to atoms, one may imply that the fc = mode of the 
photon develops a ground state, with the expectation value 



{O'sym,k—o} — 



\O'ant,k=0) — 



9sym Jsym 
^sym{k 0) 

Qant J ant 
nant{k = 0)' 



4'A — [{asym,k=a) + (aant,fc=o)] ) 

V'B — [{0'sym,k=i)) — {O'ant,k=o)] ■ 



(25) 
(26) 
(27) 
(28) 
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Eqs. (|25ll28p together with Eq. (ITO)) are the main resuhs of the present work, which form a 
self-consistent set of equations determining the photonic superfluid order parameters ^a(B) 
and the atomic coherent state order parameter for TLAs Ja(b)- Based on these results, 
the proposed two-site Hamiltonian will be used to study the QPT of light analytically and 
numerically in the following. 

In general, there are more than one solution to the self-consistent mean-field equations. 
But by calculating the corresponding free energy, one can determined the equilibrium state 
which occupies with the lowest free energy. The free energy density of our system Fs can be 
written as Fs = Ff+Fp — Em , where Ff, Fp, and Em are the free energy densities associated 
with the mean- field Hamiltonian in Eq. ([SJ, photonic mean- field Hamiltonian in Eq. (j20[) . 
and the correction term to the double counting from our mean-field decomposition in Eq. 
([3]), respectively, i.e. 



Fr, 



-Ea -Eb-^ [ln(l + e-^^^) + ln(l + e-^^«)] , 



Y jj \y 

k 



{gAJAf {gBJB? 



En 



{gA-JAY , {qbJb) 



21 



J 



A(B) 



fpA(B)gA(B) 



2E 



A{B) 



and the corresponding ground state energy density per cavity for the system is 

'{gAJAf , (gB'JB)^ 



E„ 



Ea-Ee 



(29) 
(30) 

(31) 



where V is the volume of the system. 

At zero temperature, the ground state for the mean-field Hamiltonian in Eq. ([TU| has 
the simple relation. 



(32) 



(33) 



Combine the mean-field equations in (|25II28|) and Eq. 
following matrix form. 



i'A 
i^B 



10|) . one can obtain ipA{B) in the 



(34) 



where = VL,^„^ 



n:. 



and nz^ = nzL - nz 



aiiLL — ^ugyni — ^'■ant the related photon dispersions. 
Eq. (p4)) supports an trivial solution to the mean- field equations, i.e. "0^ = "05 =0. But 
non-trivial solutions can be obtained only numerically in general for the reason that Ea(b) 
is a function of ipA{B} in Eq. ([9]). 

To give a simple picture for the QPTs of light in our system, we first consider the special 
case when the atoms A and B are identical. For this case ipA = 4'b = "0! Ea = Eb = E, 
^A ~ ^B — ^, and gA — gB = 9, can derive the photonic superfluid order parameter 
explicitly. 







g 



AEfl 



sym 



E 

Tanh{ — )il). 



(35) 



In particular, non-trivial solution, 0^0, exists at T = only when the following inequality 
is satisfied. 



g 

CI 



(36) 
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and the corresponding order parameter solution is 



9 



(37) 



which indicates a zero-temperature phase transition of superfluid to Mott-insulation phases 
at the condition 



(38) 



A corresponding finite temperature transition also exists at the critical temperature given 

by 



-^ = Tanh{—). 



(39) 



To illustrate what happens when the atoms A and B are different, we consider the limit 
case of 171^ = (or k = 0), i.e. without inter-species hopping effects. In this case, at 
T — Q the atoms A and B couple separately to the radiation fields. From Eq. (|38l) . each 
of the two species has a MI to SF phase transition at gA = \J Ail+eA and gs = \/4f7+eB, 
respectively. In such a way, there should be four possible phases in the parameter plane (g^, 
5b). These four phases correspond to (1) both of the two-species atoms are in the MI state, 
^ 4r2+e^(B); (2) only A-type atoms are in the SF state, named as the SF-A state for 
g\ < 412+6^ and g^ > 412+eB; (3) the SF-B state with the case for only i3-type atoms are 
in the SF state; and (4) all the atoms are in the SF state, i.e. the co-existence SF-AB state 

9AiB) > 4rj+e^(B). 

To give a qualitative analysis on the possible four phase states, we assume nZ^ to be 
small and perform a perturbative expansion on the limit case flZ^ = for the solutions 
of the mean- field Hamiltonian in Eq. p4p at zero temperature T = 0. By expanding the 



superfluid order parameter ipA(B) — 4'a(b) 
obtain, 



AiB)- 



to the zero-th order one can again 



'f^A{B) 







gA{B) 

4f2+ 



gA{B) 



for 
for 



3a(b) , ^ 

-mr < ^A{B) 



4n+ 



> 



^A(B) 



(40) 



and to the first order expansion in 0_ , 

/ 2 

9b(A) 

A 1 i 

Ha{B) = 



A(B) 



(0)2 



-V' 



(0) 
B{A) 



(0) 



(0)2 



(0)3 



^A{ 



for ip 



for 



(0) 
A{B) 



(41) 



with E^°lg^ = 



(0)2 
MB) 



MB) 



With above perturbative results, one can have a clear physical interpretation for the 
QPT in our proposed system. When all the atoms are in the MI phase, ■0^'' = = 0, 
both of the perturbed superfiuid order parameters S^'a and StJjb are zero as expected. The 
MI phase of the system is not modified by the perturbation in fiZ^. On the contrary, the 
properties of SF states are however strongly modified. From the second line in Eq. (PT|). 
one can easily find that 6^JJA(B) is non-zero as long as ^^b^ia) 7^ 0, which is independent of 



whether V'^(b) zero or not. In such a scenario, one of atomic species in the superfluid 
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Fig. 2: Phase diagram on the parameter plane {qa and qb) at zero temperature, i.e. for different 
atom-photon coupHng constants (in the unite of frequency). Four different phases are indicated 
in the plot. Other parameters used in the simulations are k = 0.4, k = 0.2, u — 2.7, ^ = 0.2, 
eA = 2.7, and es = 2.5. 

phase can induce a non-zero superfluid order parameter on the other species of atoms, which 
is originally in the MI state. Unlike the case for only one species of atoms, in our system 
both species of atoms develop nonzero order parameters with the superfluidity driven by 
another type of atomic species, similar to the Josephson effect in the two weakly coupled 
superconductors [T7] . 

The coupling between two atomic species through the radiation fields smears out the 
difference between the three superfluid phases mentioned above, i.e SF-A, SF-B, and the 
coexistence SF-AB phases, and turns the phase transitions between these different phases 
into crossovers. Nevertheless the qualitative properties of these three SF states are quite 
different. For example, when turning on an additional laser with the frequency ~ on our 
system, the input light resonates strongly with the S-type atoms, and the superfluidity of A- 
type atoms would be sequentially modified depending on their original state of the system. 
In an original SF-B phase, superfluidity of A-type atoms would be destroyed because of 
decoupling from B-type atoms, while an original co-existence SF-AB phase would be driven 
into the SF-A phase. 

To verify above perturbative analyses, we solve the mean-field Hamiltonian in Eq. (I34p 
by using direct numerical simulations. For the zero temperature, we solve the equations 
for different values of gA and gB by fixing the value of transition energy £a{b) and photon 
dispersion • Fig. [2] demonstrates the phase diagram of our system in the parameter 

plane gA and gs at zero temperature. As conjectured by the perturbative analysis, we 
have the phase diagram for the QPT of light in our dual-species configuration, where the 
boundary for SF-A/B to SF-AB phases is turned into crossovers. Due to the Josephson-like 
coupling effect mentioned above. The co-existence SF-AB state in Fig. [5] also occupies a 
much broaden area in the direct numerical simulations. 

For non-zero temperature, we plot in Fig. [Sja) and (b) the superfluid order parameters 
ipA and ipB as a function of the temperature, T, with different values of gA and gs, respec- 
tively. As shown in Fig. EJa), for the parameters gA = 0.1 and gs = 2.0, the system is in 
the SF-B phase at zero temperature. A clear flnite temperature insulator to superfluid tran- 
sition is found at the critical temperature Tc ~ 5. Notice that ipB > i^A but the differences 
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Fig. 3: Superfluid order parameters tpA and versus the temperature T (in the unit of Kelvin) 
with different values of (a) qa = 2.0, gs = 0.1; and (b) qa ~ 2.0, qb = 2.0. Other parameters used 
are the same as those in Fig. [2l Ifere the hopping constants are k — 0.4 and k' = 0.2. 

between tpj^ and V'b remains small throughout the whole temperature range in the SF-B 
phase, indicating the importance of Josephson coupling. For the co-existence SF-AB state 
at zero temperature, the two curves for ipA and ipB stay close to each other for the whole 
temperature range. The overall magnitudes of ipA and V'b are larger by a ratio of factor 2 
when compared with the case with gA = 0.1 and gs — 2.0. 

To examine more carefully the effect of Josephson-like coupling, we show again in Fig. 
13] the temperature dependence of ipA and with all the same parameters except changing 
the hopping constants to k' = 0.175 and k = 0.35. Comparing to the cases in Fig. [3l we 
find that the magnitudes of tJja and ips become both smaller due to the reduced Josepshon 
coupling effect, and the difference between ^Jja and 4'b becomes larger. Notice that the 
values of ^Isym and fiant are increased by decreasing k' and k (see equations (26), (27), (31), 
(32)) and the system is driven towards the MI phase. In fact we find that V'A = V'b = 
and the system is already in the MI regime for k' = 0.2 and k — 0.3 with our chosen set of 
parameters. 

With the two-site Hamiltonian by including two different atomic species, we show the 
phase diagrams for the quantum phase transitions of light by applying the Bose-Hubbard 
Hamiltonian for photons with the Jaynes-Cummings photon- atom interaction. Using a 
self-consistent mean-field approximation, we analyze the equations for the superfluid order 
parameters analytically and numerically. Four different phases are found in our proposed 
system, including an isolated Mott-insulator phase, and three superfluid phases labelled 
as SF-A, SF-B, and co-existence SF-AB states. The transitions between the different su- 
perfluid phases are found to be smeared out by the Josephson-like coupling effect between 
different types of atoms. Our results demonstrate the possibility to implement a photonic 
cavity system as a quantum simulator based on different atomic species. As the studies in 
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Fig. 4: Superfluid order parameters ipA and ipB versus the temperature T (in the unit of Kelvin). 
All the parameters used are the same as those in Fig. O but with different hopping constants, 
K = 0.35 and k = 0.175. 
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condensed matter physics, more exotic phases and richer phase diagrams for the quantum 
phase transitions of light are expected with more comphcated configurations and multiple 
atomic species in the proposed model. 

* * * 
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